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A possible spin-reorientation transition in field-induced magnetic ordering phases of the 
S = I Haldane system with large easy-plane anisotropy is proposed, using an effective 
Lagrangian formalism as well as the density matrix renormalization group method. Such 
a spin-reorientation transition is predicted in the case where the applied magnetic field is 
inclined from the easy axis of the anisotropy. We point out that this transition has a close 
connection with a variation of the order parameter even at zero temperature, although it is 
different from a quantum analog of the so-called spin-flop transition proposed for the system 
having a strong easy axis anisotropy. In connection with a novel phase observed recently 
in the Haldane system at high fields, we discuss possible implications for the field-induced 
magnetic ordering. 

KEYWORDS: anisotropic Haldane chain, field-induced phase transition, spin-reorientation 
1. Introduction 

Quantum fluctuations exhibit interesting roles in one-dimensional antiferromagnetic spin 
systems and eventually destroy the antiferromagnetic ordering even at zero temperature. As 
a consequence, the ground state turns out to be a novel singlet state having an energy gap to 
triplet states. An S = 1 spin chain with nearest neighbor exchange interactions, the so-called 
Haldane chain, is a typical example of such systems. 1 ) 

In order to understand such a ground state, many researchers are interested in the re- 
covery of magnetism, magnetic ordering, destroying the gap by a strong external magnetic 
field, for instance in Ni^HgNa^NC^ClOzi). 2 ) Honda et al. were the first to realize exper- 
imentally the field-induced long-range Neel ordering at about 4 T in the Haldane system, 
Ni(C 5 Hi4N2)2N 3 (PF 6 ), abbreviated hereafter as NDMAP. 3 " 6 ) It is noted that the critical 
behavior of the phase transition depends crucially on the direction of the magnetic field, sug- 
gesting an important role of the anisotropics of a Ni 2+ spin: it has a rather large easy-plane 
anisotropy with a weak in-plane anisotropy. 7 ' 

Following the work of Honda et al, Tsujii et al. extended the measurement of the heat 
capacity in fields up to 30 T along the crystallographic c axis and found a new phase transition 
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at 14 T in addition to the former at 4 T. 8 ) At first glance, this transition seems to be a quantum 
analog of the so-called spin-flop transition, which occurs in the case where the ground state 
is the Neel state. 9 ) We, however, realize immediately that it is not the case because NDMAP 
has the disordered ground state and the field-induced ordered moment is perpendicular to it. 

Now, we remember the tilting of the principal axes of the anisotropy in NDMAP: the 
local principal z axis does not coincide with the crystallographic c axis but tilts from the c 
axis by 16 degrees in the a-c plane. Furthermore, we note that there are two types of chains, 
each of which belongs to one of the sublattices of the body-centered-cubic lattice and has 
the tilt direction just opposite to another. Thus, the antiferromagnetic ordering of the two 
sublattices are magnetically decoupled at the mean-field level due to the geometric frustration. 
It is worth noticing that such an interesting situation has skillfully been confirmed recently 
by the neutron elastic scattering: a half-ordered state has been observed in NDMAP between 
the two critical fields for each sublattice, adjusting the field direction to the z axis of one of 
the two kinds of sublattices. 10 ) 

In order to understand what happens in NDMAP at high fields, we will consider the 
following Hamiltonian describing the Haldane system with the planar anisotropy as well as a 
weak in-plane anisotropy under the external magnetic field, 

H/J = £ [S< ■ S l+1 + d{Sf ) 2 - e{(Sf f - (Sf) 2 } -h-Si], (1.1) 

i 

where S is the S = 1 spin operator and the Cartesian coordinates (x, y, z) are referred to 
the spin anisotropy axes. The parameters, J, d, e and h denote, respectively, the exchange 
constant, which defines the unit of energy throughout this paper, the easy-plane anisotropy 
constant d(= 0.25), the in-plane anisotropy constant e(= 0.01) and the external magnetic 
field. u ~ 13 ) Note that the anisotropies d and e make the x axis easy while the z axis hard. 
Considering that the anisotropy axis in NDMAP tilts in the a-c plane, 10 ) we assume h has 
the x and the z components: h = /i(sin6>, 0, cos #). A similar system has been studied by the 
perturbative calculation. 14 ) The aim of this paper is to explore the role of quantum fluctu- 
ations in the field-induced ordering phase from the unified point of view. Especially, we will 
be concerned with the nature of various ordering phases and phase transitions between them, 
focusing our attention to the field direction. We are interested in the field region above the 
critical field h c , at which the Haldane gap is closed. To this end, we adopt the density ma- 
trix renormalization group (DMRG) method 15 ) as well as the phenomenological field theory 
(PFT); 16 ) the former is used for ab initio simulations for a single spin chain while the latter 
is suitable for unified understanding of the phenomena. The former is also used for the con- 
firmation of the results of PFT. With the aid of the information on the single Haldane chain, 
we analyze the phase diagram, using the mean-field approximation for weak inter-chain inter- 
actions. It has been confirmed that such an approximation reproduces the phase diagrams of 
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NDMAP quite well. 4 ) It is worth noticing that the important information on the criticality of 
the single spin chain can be drawn from the analysis of the phase diagram. 4 ' 17 ) 

The present paper is organized as follows. In the next section, we describe the methods 
used to obtain the field-induced ordering states in the Haldane system and their excited states. 
In §3, we give theoretical results and consider the phase diagram and other thermodynamic 
quantities. In the last section, we discuss possible implications and relevance of the present 
results in connection with the experimental results in NDMAP and conclude our discussion 
with some remarks. 

2. Methods 

In this paper, we perform simulations for the Haldane phase and ordering phases in the 
anisotropic Haldane chain, using (1) the ab initio numerical calculation by the DMRG method, 
and (2) the Landau-Ginsburg type theory of PFT, within the mean-field approximation. As 
will be seen soon later, the direction of the applied magnetic field with respect to the anisotropy 
axes plays important roles in the physics of the field-induced phase transition of the anisotropic 
Haldane chain. 

In the ordering phase, we define the sublattice magnetization 

L = (S>2n — S^n+l) , (2.1) 

and characterize possible ordering phases in the h-0 phase diagram in addition to the Haldane 
phase: if L is in the x-z plane, we call it the xz-phase and if L is in the y direction, we call 
it the y-phase. 

2.1 Density matrix renormalization group method 

The DMRG method 15 ^ has been developed in order to reveal properties of one-dimensional 
strongly correlated systems. It is known that this method provides the ground state energy 
and low-lying excited state energies of spin-gaped systems with very high accuracy. Then, we 
use the DMRG method to calculate the physical quantities of the Haldane system. Further, 
we extend it so that it is applicable to the ordering states, and discuss the phase transition 
between the xz-phase and the y-phase. 

The anisotropic Haldane chain described by the Hamiltonian (1.1) is studied by the 
infinite-system DMRG algorithm. 15 ) In this algorithm, the whole system is divided into two 
blocks, system and environment, as is schematically shown in Fig. 1. These blocks are elon- 
gated iteratively keeping the Hilbert space up to a certain size m. The bases of the system 
block of length I + 1 are chosen from eigenstates with the largest eigenvalue of the reduced 
density matrix: 

p s = Tr E |V> (VI , (2.2) 

where TVe means tracing out the degrees of freedom of environment side and ip is the ground 
state wave function of the superblock of length 21 + 2, which is built from the system and 
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Fig. 1. A schematic diagram of the superblock. The superblock of length 21 + 2 is constructed from 
the system and the environment blocks of length I and two single sites. Blocks of length l + l, Sj+i 
and E; + i, are obtained from the superblock of length 21 + 2. 



environment blocks of length I and two single sites between them. Usually, we can form the 
environment block by reflecting the system block. 

The calculation has been performed under the open boundary condition (OBC). As is well 
known, since the Haldane system with OBC has two effective 5 = 1/2 degrees of freedom at 
each end of the chain. 18 ) Then, we add a real 5 = 1/2 spin at each end in order to remove 
the ground state degeneracy. 

In order to study the ordering phase, we need to extend DMRG so that it is able to treat 
a phase in which the Z2 symmetry is broken. To this end, we first apply a weak magnetic field 
to an edge spin of the superblock to break the Z2 symmetry. As a consequence, we choose 
the ground state \ip) of the superblock uniquely. Then, we construct the system block and the 
environment block, using eq. (2.2) individually for each block. Repeating these procedures, 
we finally obtain the ordering state. 

We compute the energy gap A, which is obtained by subtracting the ground state energy 
of the superblock from the lowest excited state energy, and the sublattice magnetization L in 
the ordering phase, which is evaluated at two spins located on the center of the superblock, as 
a function of the field h at various values of the field angle 9. As the infinite-system DMRG 
is often trapped into metastable states, the calculation has been done starting both from 
the xz-phase to the y-phase and from the y-phase to the xz-phase at the vicinity of the 
transition point. These calculations show the hysteresis, which is reasonable for the first order 
phase transition. It is, however, noted that the hysteresis region of the field observed in our 
calculations is wider than usual, on which we will discuss later. 
2.2 Phenomenological field theory 

Using a Landau-Ginsburg type field theory, Zheludev et a/. 19 ) have studied the transition 
between the Haldane phase and the ordering phase in NDMAP. They succeeded in explaining 
some aspects of the transition. We follow their method but extend their Lagrangian, so that it 
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can describe the phase transition between the ordering phases: we need more care about the 
anisotropy in higher order terms. We start with the S = 1/2 alternating chain Hamiltonian, 19 ) 



n 



E 



(2.3) 



where J a , J' a and h mean, respectively, two kinds of nearest-neighbor anisotropic exchange 
constants and the external magnetic field. In comparison with the Hamiltonian they have 
used, this Hamiltonian is modified in the J' term where the anisotropy is introduced. The 
Lagrangian obtained by the similar procedures as described in ref. 19, is written as 

C = E f($&) 2 " v i - ( h x 0)i ~ U * (0) " ^ (0, $0) , (2.4) 

i=x,y,z ' 

where v and denote, respectively, the staggered order parameter, which is related to 
the sublattice magnetization, the characteristic velocity and the external magnetic field. The 
quadratic and quartic parts of the potential are given by 



rn 



1 



rrv. 



(h x 4>)\ 



u 4 (0, d t <j>) = E Ml<t> 2 + E Kh&^f] + E /•;/•)• 



(2.5) 
(2.6) 
(2.7) 



F = -d t cj) + h x <f). 

We note that the modification appear only in U4, which becomes more anisotropic. If we set 
Aj = A, X± t ij = Ai, A2,jj = A2, we obtain the quartic potential obtained by Zheludev et al. 19 ^ 
Note that \2,ij is a symmetric tensor, A2,ij = A2ji, since 8^(8^1/4) = 8^(8^1/4). 

The transition from the Haldane phase to the ordering phase is suggested by a new static 
non-zero solution, where the static part of the potential has a minimum. This transition 
is characterized by the change of the sign of the second derivatives of U<i- Two kinds of 
transitions, namely transitions from the Haldane phase to the xz-ph&se and to the y-phase, 
are possible depending on 9. The critical field h c as a function of 6 is obtained as 

h c = mm(hi xz \h^y (2.8) 



hi xz) = 



m x m y m z 



m x sin 6 + m z cos 2 9 ' 



m x m y m z 



(2.9) 



(2.10) 



y m x sin 9 + m z cos z 9 

where lif^ and are the critical fields for two transitions. Note that lif^ and are 
independent from the modification of our Lagrangian since they depend only on Ui. 
The static solution (f>^ can be obtained from the following simultaneous equation: 

E ttij<pj + E ^U,mn<Pl<Pm4>n = 0, (2.11) 



Imn 
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where the matrices f2, A are defined as 

Qjj = mjSjj — y~] eun e jl'n ~ 1 i (2.12) 

ll'n 

A-ij,mn = ^ij,mn ~l~ ^mn,ijj (2-13) 
^ij,mn = \$ij$mn + ^ ] ^l,ik^ij^kk' + ^2,ij^ik^jk'^) ^kln^k'l'm~ ~ • 

(2.14) 

Here e^-fc means the Levi-Civita symbol. 

The excitation energies are obtained by the equation of motion, where the small fluctua- 
tions around the static solution are treated up to quadratic order. There are three excitation 
modes, of which energies are obtained as three real roots of the following secular equation: 



det (M - uj 2 G - iwC) = 0, (2.15) 



where M, G, C are defined by 



Mij = (Uj + ^-Sij + ( A «3,™ + A im,jn + Kn,mi) <f>$ 'ffl , (2-16) 
TYli 



Gij 



m ■ ^ — ^ 



ra ' i 




4°v; o) 

Sij - \2 tij -^-, (2.17) 



miirij 

Cij = Rij — Rji, (2-18) 

Rij = ^r- E + E {^l,k'i^ilkSjk' + ^l,k' j^ilj^kk 1 ) ^^(fi 

I 

(0) 



mi I kk'l 



nii *-f V m 3 - m k J J 



(2.19) 



The phenomenological parameters, m, m, A, Ai, and A2, are chosen so that PFT reproduces 
the magnetic field dependence of the excitation energies observed by the neutron inelastic 
scattering (NIS) and by the electron spin resonance (ESR) experiments under the magnetic 
fields along the a and the c axes. Note that the case in which the magnetic field along the a 
(c) axis corresponds to the case 9 = 0.67T (9 = O.lir) for a spin chain and 9 = OAtt (9 = — Q.Itt) 
for another chain in NDMAP. We determine six parameters, m and m, from the data of the 
energy gaps and of the critical fields in the Haldane phase. Remaining parameters, A, Ai and 
A2, are chosen so that PFT reproduces the excitation energies in the ordering phase (h > h c ). 
Although these parameters can not be determined uniquely from NIS and ESR results, we 
chose a following set of the parameters with the aid of the result by DMRG, which will be 
shown soon later: 

m x = 0.19, m v = 0.29, m z = 1.7 

' (2.20) 
m x = 0.15, rhy = 0.15, rh z = 0.32 
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Fig. 2. The lowest three excitation energies as a function of the magnetic field along (i) the c axis 
and (ii) the a axis. The circles and the squares denote the results observed by ESR 20 ^ and by 
NIS 21 ) with the wave vector q = n, respectively. The result calculated by PFT is shown by the 
solid lines. The vertical dashed lines denote the critical field h c , which is the boundary between 
the Haldane phase and the ordering phase. 



A = (1.4 1.0 1.0) , 



Ai 



/ 0.37 0.37 0.37 \ 

0.37 0.37 0.25 
v 0.37 0.48 0.37 j 



( 0.30 0.30 0.49 \ 

0.30 0.30 0.44 
v 0.49 0.44 0.30 y 



(2.21) 



In Fig. 2, we show the magnetic field dependence of the excitation energies calculated by PFT, 
which are compared with the experimental results obtained by ESR 20 ) and by NIS. 21 ) Our 
PFT can reproduce quite well the experimental results. 

3. Results 

3.1 Energy Gap and Sublattice Magnetization 

The energy gap A and the sublattice magnetization L in each phase are calculated as 
a function of the magnitude h and the angle from the z axis 9 of the magnetic field, using 
DMRG as well as PFT. The result calculated by DMRG as a function of h at the fixed angles, 
9 = and 9 = tt/2, is shown in Fig. 3. In both cases, the staggered ordering state appears 
above h c . The direction of the sublattice magnetization, however, is different, depending on 
9: L is parallel to the x axis for 9 = 0, while it is parallel to the y axis for 9 = tt/2. The latter 
case corresponds to the situation of the spin-flop transition in the classical spin system but, 
in our quantum case, it is not difficult to realize that there is no quantum analog since the 
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Fig. 3. The energy gap A and the sublattice magnetization L as a function of h at (i) = and (ii) 
9 = tt/2 calculated by DMRG. The calculation has been done keeping up to 45 states (m = 45). 
The Haldane gap closes in the case (i) at h c = 0.26, and in the case (ii) at h c — 0.45, and a gap 
open again in the ordering phase. The field-induced ordering phase is the x-phase for 8 = 0, and 
the y-phase for 9 — tt/2. 



induced staggered magnetization is already perpendicular to the field. 

In Fig. 4, we show the result calculated by PFT. Because our Lagrangian describes be- 
havior of the system only in the small values of we show the result up to h = 1. The 
order parameter is parallel to the x axis for = and to the y axis for 9 = tt/2, which is 
qualitatively the same as the sublattice magnetization obtained by DMRG in Fig. 3. 

Next, let us consider the cases of intermediate angles. We show A and L (4>) calculated by 
DMRG and by PFT for = O.lir, whose situation corresponds to the case where the magnetic 
field is applied along the c axis in NDMAP. The results obtained by DMRG and by PFT are 
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(i) = 

0.2 | . — — ■ . ,— — , 1 | . . . . 1 




Fig. 4. The energy gap A and the order parameter cf> as a function of h at (i) 6 = and (ii) 9 = it/2 
calculated by PFT. The Haldane gap closes in the case (i) at h c — 0.17, and in the case (ii) at 
h c = 0.30. Then, the field-induced ordering phase occurs. The orientation of the order parameter 
corresponds to the sublattice magnetization in DMRG; the x-phase for 9 = and the y-phase for 
9 = tt/2 

qualitatively consistent: with increasing the magnetic field, the Haldane gap becomes smaller 
and closes at h c , where the quantum phase transition occurs from the Haldane phase to the 
xz-phase similar to the case 9 = 0. It is surprising to find a successive phase transition: the 
xz-phase further changes to the y-phase at the higher critical field, h sr . Since this transition 
is first order, the xz-phase coexists with the y-phase in the field region h ST i < h < h ST 2, where 
h S ri (h ST 2) is the field at which the y-phase (the xz-phase) becomes unstable. The latter phase 
transition between two kinds of the ordering phases can be said to be "spin-reorientation 
transition" . 
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Fig. 5. The energy gap A and the sublattice magnetization L (cf>) as a function of h at 6 = O.ln 
calculated by DMRG and by PFT. Two successive phase transitions can be seen in these figures: 
the one is the order-disorder phase transition at h c , and the other is the first order phase transition 
from the xz-phase to the j/-phase: the spin-reorientation transition. The xz-pha.se and the j/-phase 
become unstable, respectively, at h ST 2 and at h ST \. 

Next, we show A and L (0) as a function of 9, in order to investigate the spin-reorientation 
transition. In Fig. 6, we see that, with increasing the field angle 9, the xz-phase becomes 
unstable against the y-phase. At the intermediate angle, # sr i < 6 < 9 ST 2, both phases coexist, 
which is a typical signature of the first order phase transitions. 
3.2 Phase Diagram in the h-9 plane 

We calculate the energy gap and the sublattice magnetization with respect to h and 9 by 
DMRG as well as by PFT. From these results, we extract the phase diagram in the h-9 plane 
shown in Fig. 7. The phase diagrams obtained by DMRG and by PFT do not coincide with 
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Fig. 6. The energy gap A and the sublattice magnetization L (<fi) as a function of the field angle 
at h — 0.5 by DMRG and by PFT. The xz-phase and the y-phase become unstable, respectively, 
at # sr 2 and at 9 sr i. 



each other but are qualitatively similar: These phase diagrams consist of the Haldane phase, 
the xz-phase and the y-phase. It is interesting to note that when the applied magnetic field 
tilts from the z axis, the spin-reorientation transition can occur at the higher critical field h sr 
following the order-disorder phase transition at the lower critical field h c . 

4. Discussions 

We have studied the field-induced phase transition in the anisotropic S = 1 Haldane chain 
and found the new phase transition at the higher critical field than that of the order-disorder 
phase transition by means of DMRG and of PFT. Under the external magnetic field in the 
x-z plane, depending on the field angle 9, two kinds of the ordering phases are induced at 
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(a) DMRG (b) PFT 
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Fig. 7. The phase diagram in the h-9 plane obtained (a) by DMRG and (b) by PFT. In (a), the 
closed and the open circles represent the critical fields, respectively, from the Haldanc phase to 
the xz-phase and from the Haldane phase to the y-phase, while the closed and the open diamonds 
denote, respectively, the points, where the xz-phase and the j/-phase become unstable. The lines 
are guides for the eye. In (b), the solid and the dashed curves correspond, respectively, to the 
closed and the open circles in (a), while the dotted and the chain curves correspond, respectively, 
to the solid and the open diamonds. 



h c , at which the Haldane gap closes: when the magnetic field is applied along near the z axis 
the xz-phase is induced, while it is along near the x axis the y-phase is induced. The spin- 
reorientation transition from the former to the latter further occur in a certain range of the 
angle 9. In the vicinity of the spin-reorientation transition field, the xz-phase and the y-phase 
coexist together and the magnetization process shows hysteresis, which is typical for the first 
order phase transition. 

Although our result of DMRG keeping the states m = 45 in the cycle at fields near 
the higher critical field does not converge so well with respect to m, it does in fields near 
the lower critical field. Then, for a quantitative comparison of the high field transition with 
the experimental result, we need more sophisticated calculations, although at present we 
can not extend m in our calculation because of the computational limitation. On the other 
hand, PFT is flexible involving many parameters but it is phenomenological. Nevertheless, 
we may expect that these methods provide a qualitatively correct description, considering the 
reasonable results obtained by both methods. 

Considering these limitations, let us remind the experimental finding in NDMAP, a new 
phase transition at the high field. The novel phase transition of the first order at high fields 
predicted by our theory is favorably compared with the phase transition observed, since in 
the experiment the magnetic field is applied along the c axis {\9/tt\ = 0.1 for both sublattices) 
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in NDMAP is located in the field region where the successive phase transition occurs. We 
mention that this is a consequence of the tilting of the anisotropy axes from those of the 
crystal. It is interesting to confirm experimentally the fact predicted: the critical field h sr is 
very sensitive to the angle 6 near 6/ir ~ 0.1 as is seen in Fig. 7. Also NMR experiments to study 
the direction of the order parameter in each phase are desirable. Another point we should note 
for the application of our theory to NDMAP is a choice of the anisotropy parameters d = 0.25 
and e = 0.01. These values, especially the latter, have not been determined so accurately yet, 
but it dominates the critical field of the spin-reorientation transition. Then, in order to study 
this transition more quantitatively it needs to adjust these parameter values more carefully 
considering other experimental results. 

In conclusion, with the aid of the results of DMRG as well as of PFT, we can draw the 
physics in the field-induced quantum ordering phases and that in the phase transition between 
them, which gives rise to the scenario to explain the experimental finding in NDMAP, the 
new phase transition at the high field. 
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